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We review the status of the calculation of next-to-next-to-leading order corrections to large angle Bhabha 
scattering in pure QED. After discussing the electron- loop and photonic corrections, we focus on the recently 
calculated two-loop virtual corrections involving a heavy-flavor fermion loop. We conclude by assessing the 
numerical impact of these corrections on the Bhabha scattering cross section at colliders operating at a center of 
mass energy of about 1-GeV. 



1. INTRODUCTION 

The scattering process e~ e + — * e~ e + is named 
after the author that first calculated its cross sec- 
tion in 1936 pp. Nowadays, Bhabha scattering 
is of crucial importance in the study of elemen- 
tary particle phenomenology, since it is the pro- 
cess chosen for the luminosity determination at all 
e + e~ colliders. In particular, at colliders operat- 
ing at a center of mass energy y/s ~ 1 — 10 GeV, 
the luminosity is measured by observing Bhabha 
events at large scattering angles. In this kine- 
matic region, the Bhabha scattering cross sec- 
tion is large and QED dominated: these two 
properties allow precise experimental measure- 
ments and an accurate theoretical evaluation of 
the cross section in perturbation theory. Any re- 
alistic calculation of the theoretical Bhabha scat- 
tering cross section must take into account the 
detector geometry and experimental cuts. For 
this purpose a number of sophisticated MC gen- 
erators were developed [2|3|4|5|6] , In particu- 
lar, at KLOE/DA<I>NE, the experimental accu- 
racy reached on the measurement of Bhabha scat- 
tering events, joined with the theoretical accuracy 
of the MC generator BABAYAGA (3I4I7| . allowed to 
determine the luminosity with an error of 0.1% 
[8j. Since the theoretical error on the Bhabha 
scattering differential cross section directly af- 
fects the luminosity measurement precision, in 
recent years a significant effort was devoted to 



the calculation of perturbative corrections to this 
process. The next-to-leading order corrections 
are well known even within the full Standard 
Model [9]. At the next-to-next-to- leading order 
(NNLO), for what concerns the electroweak cor- 
rections, only the logarithmically enhanced terms 
were calculated [10]. In pure QED, instead, the 
situation is sensibly different. The first diagram- 
matic calculation of the two-loop QED virtual 
corrections to Bhabha scattering can be found in 
|llj . However, this result was obtained by em- 
ploying dimensional regularization (DR) to reg- 
ulate both soft and collinear divergencies, while 
all the available MC generators employ the elec- 
tron mass m e as a collinear regulator. Today, 
the complete set of NNLO corrections to Bhabha 
scattering in pure QED have been evaluated using 
m e as a collinear regulator. The two- loop Feyn- 
man diagrams involved in the calculation can be 
divided in three gauge independent sets: i) di- 
agrams without fermion loops ("photonic" dia- 
grams), ii) diagrams involving a closed electron 
loop, and Hi) diagrams involving a closed loop of 
a fermion heavier than the electron. In the fol- 
lowing three sections we consider the calculation 
of each set of corrections separately. In Section [3 
we conclude by briefly discussing the numerical 
impact of these corrections on the Bhabha scat- 
tering cross section at y/s ~ 1 GeV. 
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2. PHOTONIC CORRECTIONS 

All the NNLO photonic corrections, with the 
exception of the ones arising from the two-loop 
photonic boxes, can be obtained in a straightfor- 
ward way, retaining the full dependence on m e 
[12] , by using a technique that is by now stan- 
dard in multi-loop calculations. This technique 
is based on the Laporta algorithm [13] for the 
reduction of the Feynman diagrams to the Mas- 
ter Integrals (Mis), and then on the differential 
equation method [XT for their analytic evalua- 
tion. However, a calculation of the two-loop pho- 
tonic boxes retaining the full dependence on m e 
is beyond the reach of this technique. This is 
due to the fact that the number of Master In- 
tegrals belonging to the same topology after the 
reduction process is larg^]. This means that one 
should be able to solve analytically large systems 
of first-order ordinary differential equations; this 
is in general not possible. Alternatively, one could 
use the Mellin-Barnes techniques [15] to calculate 
the Mis one by on^l The calculation, however, 
is very complicated and, moreover, from a phe- 
nomcnological point of view, to keep the full de- 
pendence on the electron mass is not important. 
In fact, the physical problem exhibits a well de- 
fined mass hierarchy. The mass of the electron 
is always very small compared to the other kine- 
matic invariants and it can be safely neglected 
everywhere, with the exception of the terms in 
which it acts as a collinear regulator. It can be 
shown that the collinear structure of the NNLO 
corrections is the following 
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where L e = In (s/m^) and where the coefficients 
<jW are functions of the scattering angle 9 and, 
in general, of the mass of the heavy fermions in- 
volved in the virtual corrections. For any practi- 
cal purpose, the approximation given by Eq. ([T]) 
is sufficient for a phcnomcnological description of 
the procesfl In the photonic case, the coeffi- 

2 Some topologies have up to 6 Mis. 
3 Partial results can be found in |16I17| 
4 It can be shown that the terms suppressed by a posi- 
tive power of m^/s do not play any phenomenological role 



cients of the squared and single collinear loga- 
rithm in Eq. ([T]) were obtained in [18119] , The 
precision required for luminosity measurements 
at e + e~ colliders demanded the calculation of the 
non-logarithmic coefficient, that was obtained in 
[20] through the infrared matching to the mass- 
less approximation. The technique of [20] allowed 
to reconstruct the photonic differential cross sec- 
tion in the s 3> ^ limit from the calculation 
in |11] , where m e was set to zero from the start. 
The method employed in [20] involves a change of 
regularization scheme for the collinear divergen- 
cies originating from a vanishing electron mass. 
A method based on a similar principle was subse- 
quently developed in |21I22| ; the authors of [22] 
confirmed the result of [20] for the NNLO pho- 
tonic corrections to the Bhabha scattering differ- 
ential cross section. 

3. ELECTRON LOOP CORRECTIONS 

The NNLO electron loop corrections arise from 
the interference of two-loop Feynman diagrams 
with the tree-level amplitude as well as from the 
interference of one-loop diagrams, as long as one 
of the diagrams contributing to each term in- 
volves a closed electron loop. This set of cor- 
rections presents a single two-loop box topology, 
and it is therefore technically less challenging to 
evaluate with respect to the photonic correction 
set. The calculation of the electron loop correc- 
tions was completed a few years ago [23J ; the final 
result retains the full dependence of the differen- 
tial cross section on the electron mass m e . The 
Master Integrals involved in the calculation were 
identified by means of the Laporta algorithm and 
evaluated with the differential equation method. 
As expected, after UV renormalization the dif- 
ferential cross section presented only residual IR 
poles that were removed by adding the contribu- 
tion of the soft photon emission diagrams. The 
resulting NNLO differential cross section could 
be conveniently written in terms of 1- and 2- 
dimcnsional Harmonic Polylogarithms (HPLs) of 

already at very low cm. energies, yfs ~ 10 MeV. More- 
over, the terms rn^/t (or m 2 /u) become important in the 
very forward (backward) region, unreachable for the ex- 
perimental set ups. 
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maximum weight three [53] . Expanding the cross 
section in the limit s,\t\ ^> m^, the ratio of the 
NNLO corrections to the Born cross section can 
be written as in Eq. fTJ): 



dcr (Bo t „) 




The explicit expression of all the coefficients 
<$( EL >0 5 obtained by expanding the results of [53] 
was confirmed by two different groups |22l25j . It 
is easy to check that the cubed collinear logarithm 
in Eq. @ cancels against the corresponding term 
arising from the soft pair production graphs |26] . 

4. HEAVY FLAVOR CORRECTIONS 

Finally, we consider the corrections originat- 
ing from two-loop Feynman diagrams involving a 
heavy flavor fermion loop0. Since this set of cor- 
rections involves one more mass scale with respect 
to the corrections analyzed in the previous sec- 
tions, a direct diagrammatic calculation is in prin- 
ciple a more challenging task. Recently, Becher 
and Melnikov applied their technique based on 
SCET to Bhabha scattering and obtained the 
heavy flavor NNLO corrections in the limit in 
which s, |i|,|u| 3> rnj ^S> m^, where mj is the 
mass of the heavy fermion running in the loop 
22J. Their result was very soon confirmed by 
means of a method based on the asymptotic ex- 
pansion of Mellin Barnes representation of the 
Master Integrals involved in the calculation [25 . 
However, the results obtained in the approxima- 
tion s, 1 1 1, \u | 3> mj ^> m 2 e cannot be applied to 
the case in which the y/s < rrif (as in the case 
of a tau loop at y/s ~ 1 GeV), and they apply 
only to a relatively narrow angular region per- 
pendicular to the beam direction when yfs is not 
very much larger than m/ (as in the case of top 
quark loops at ILC). It was therefore necessary to 
calculate the heavy flavor corrections to Bhabha 
scattering assuming only that s,\t\, \u\, m 2 * 3> m% 
The technical problem can be simplified by con- 
sidering carefully, once more, the structure of the 

5 In this context, and at \fs ~ 1 GeV, by heavy flavor 
fermion we mean muon and tau leptons, as well as b- and 
c-quarks. Top quarks completely decouple at intermediate 
energies. 



collinear singularities of this set of corrections. 
The ratio of the NNLO heavy flavor corrections 
to the Born cross section is given by 




It is possible to prove that, in a physical gauge, all 
the collinear singularities factorize and can be ab- 
sorbed in the external field renormalization [27 . 
This observation has two consequences in the case 
at hand. The first one is that box diagrams are 
free of collinear divergencies in a physical gauge; 
since the sum of all boxes forms a gauge indepen- 
dent block, it can be concluded that the sum of 
all box diagrams is free of collinear divergencies in 
any gauge. The second consequence is that the 
single collinear logarithm in Eq. ([3]) arises from 
vertex corrections only. Moreover, if one chooses 
on-shell UV renormalization conditions, the irre- 
ducible two loop vertex graphs are free of collinear 
singularities. Therefore, among all the two-loop 
diagrams contributing to the NNLO heavy fla- 
vor corrections to Bhabha scattering, only the re- 
ducible vertex corrections are logarithmically di- 
vergent in the m e —* limhp, and they are very 
easy to calculate even if they depend on two dif- 
ferent masses. By exploiting these two facts, we 
were recently able to obtain the NNLO heavy fla- 
vor corrections to the Bhabha scattering differ- 
ential cross section [28129] . assuming only that 
s, \t\,\u\,mj 3> vnr e . In particular, in obtaining 
the analytic expression for the NNLO cross sec- 
tion, we worked in the Feynman gauge, setting 
m e = from the start in all the diagrams with 
the exception of the reducible ones and of the 
interference of one-loop graphs. This procedure 
allowed us to effectively eliminate a mass scale 
from the two-loop boxes, so that the evaluation 
could be carried out with the techniques already 
employed in the diagrammatic calculation of the 
electron loop correction^. We want to stress that 
in this approach individual box diagrams are sin- 
gular in the m e — > limit and the collinear singu- 
larities appear as additional poles in the dimen- 

6 Additional collinear logarithms arise also from the inter- 
ference of one-loop vertex and self-energy diagrams. 
7 The necessary Mis can be found in |29!30l 
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sional regulator e; however it is easy to prove that 
such divergencies cancel in the sum of all the box 
diagrams. By expanding the analytic results of 
28 2!j] it was possible to check the heavy flavor 
cross section in the s, \t\, \u\ 3>m^3> rv? & limit, 
which was previously known [22 25 . At inter- 
mediate energy colliders like DA<I>NE, the exact 
dependence on m/ of the results of [28 29 allows 
to account for the contribution of muons, taus, b- 
and c-quark loops to the Bhabha scattering cross 
section. In the case in which the heavy flavor 
fcrmion is a quark, it was straightforward to mod- 
ify the calculation of the two-loop self-energy dia- 
grams to obtain the mixed QED-QCD corrections 
to Bhabha scattering [35] . An alternative numeri- 
cal approach to the calculation of the heavy flavor 
corrections to Bhabha scattering, based on dis- 
persion relations, was pursued in [31]. The latter 
method allows also to evaluate the contribution 
of the light quarks vacuum polarization to the 
Bhabha scattering cross section; to this purpose 
one has to convolute the kernel functions with the 
data concerning the cross section of the process 
e + e~ — > hadrons. 

5. CONCLUSIONS 

The numerical impact of the photonic and elec- 
tron loop QED corrections to the Bhabha scatter- 
ing cross section at flavor factories was carefully 
examined in [41 7] . The authors were able to show 
that the event generator BABAYAGA, based of 
the matching of exact NLO corrections with the 
Parton Shower algorithm, has a theoretical accu- 
racy of the order of 0.1%. A similar analysis of the 
heavy flavor NNLO corrections is not yet avail- 
able. However, it is possible to evaluate numer- 
ically the NNLO heavy flavor corrections to the 
Bhabha scattering cross section. In Table Q] we 
show the results of such an evaluation for y/s = 1 
GeV and for a scattering angle 9 in the range 
50° < 6 < 130° (see [25] for details). It is possi- 
ble to observe that the muon corrections are an 
order of magnitude larger of the corrections in- 
volving heavier fermions while they are one order 
of magnitude smaller than the electron loop cor- 
rections and they reach 1/2 permille of the Born 
cross section at large scattering angles. 



In conclusion, the calculation of the two-loop 
corrections to Bhabha scattering in QED is now 
complete. The calculation of the heavy fermion 
NNLO corrections allowed to remove the last 
piece of pure theoretical uncertainty in the lumi- 
nosity determination at low-energy accelerators. 
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Table 1 

The second-order electron, muon, c-quark, r-lepton, and &-quark QED contributions to the differential 
cross section of Bhabha scattering at ^/s = 1 GeV in units of 10 -4 of the Born cross section. 
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